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The integral and local compaction occurring during sintering of ceramic materials are analyzed. A model de¬ 
scribing the structure of the material during the sintering process taking account of clusterization of the com¬ 
ponents of the particle material is proposed. Relations and an algorithm for calculating the micro- and 
macroporosity space (coordination number, specific surface, and average pore size) are obtained. Calculations 
performed with this model showed that the macropore size increases sharply at the beginning of the sintering 
process. 
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Sintering of ceramic and porous ceramic materials deter¬ 
mines their structure and properties, i.e., article quality. 

The material is not uniform during sintering because of 
clusterization, i.e., local compaction, in certain zones. The 
interconnection of these zones breaks down partially or com¬ 
pletely. As a result the macroscopic volumetric changes lag 
behind the changes in the local parts of the system [1 - 6]. 
Clusterization is characteristic mainly for the initial stage of 
sintering, and its intensity is determined by the degree of 
nonuniformity of the initial porous structure [5, 6] and the ki¬ 
netic process [3-5]. It manifests as an increase in the aver¬ 
age pore size, which decreases the shrinkage of the material 
[5]. In addition, the work in [7] should be taken into account 
indirectly. According to [7] the change in the structure of the 
solid phase of a material and the porous structure mainly pro¬ 
ceeds in different directions (as the porosity decreases the co¬ 
ordination number Z increases in the solid phase of a mate¬ 
rial and decreases in the porous structure). In the process a 
biporous structure forms during the sintering process 
[4, 8, 9], 

The sintering process for powders has been analyzed in 
[5], where the magnitudes of the local and macroscopic 
compactions were determined from the shrinkage of a porous 
body as a whole and the average pore size. In [5] it was as¬ 
sumed that at the first (hypothetical) stage the initial number 
of pores is conserved and the shrinkage AV/V m (where ARis 
the volume change and V m the initial volume of the material) 
is due to a decrease in the pore size. At this stage uniform 
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compaction (sintering) promotes a uniform reduction of the 
diameter of all pores. 

At the second stage the total pore volume is conserved 
(there is no additional macroscopic shrinkage), and the aver¬ 
age pore radius r increases from the initial value r m to the ex¬ 
perimentally determined value with the corresponding de¬ 
crease in the number of pores. This stage is due to additional 
redistribution of the material with a reduction of the free sur¬ 
face but no macroscopic shrinkage. 

However, in the method the structural model of the mate¬ 
rial in the form of uniformly distributed isolated spherical 
pores of the same size is idealized. Such a model does not re¬ 
flect the real structure of a porous material and the character 
of the change in the structure, which is an impediment to its 
use. 

A more accurate model in the form of a three-dimen¬ 
sional capillary network with average radius r taking account 
of the cluster porosity is more expedient. 

For this three-dimensional model it is possible to calcu¬ 
late the pre-sintering shrinkage of the material from the ini¬ 
tial porosity of the molded green body n m to the post¬ 
sintering porosity IT and to determine the average post-sin¬ 
tering size r of the pores (capillaries) according to the known 
initial value r m . 

A unit volume of the initial system consists of a solid 
phase V s and pores V = nnL m r* , where L m and n are, re¬ 
spectively, the average length and number of individual 
pores (capillaries) in this volume before sintering. 

The porosity decreases during the sintering process. The 
relative shrinkage of the material AV/V m is given by 
AV/V m = (IT - rr„, )/(l - n). We note that the pore sizes and 
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the shrinkage of the material are interrelated nonlinearly 
owing to the clusterization process. 

We shall examine the calculation of the change in the 
structural parameters of powder material during sintering for 
a monodisperse material comprised of particles with diame¬ 
ter D. The average structural parameters (porosity If, average 
coordination number Z (for the particles of the material), spe¬ 
cific surface S a of the material (free surface to volume ratio 
of the porous material), size d of the interparticle contacts 
and the average pore radius r) are determined according to 
the porosity of the material II, the normal local deformation 
h/D of the particles in the contact zone and the relations pre¬ 
sented in [6]: 

I — II = (1 -2/Z) 2 /[8(1.077 — 1/Z-Z U6 )(0.5 -h/D) 3 ]; (1) 

S 0 =6{\-U){\-Z[h 0 /D- d/D(h 0 /D - h/D)]}/D ; (2) 

r ~ 2T1/S 0 , (3) 

where h = h o - 4h 2 (l.5D - h o )/(3d 2 ), h o /D = [1 - (1 - 
d 2 /D 2 ) 0 5 ]/2; h is the height of a spherical segment per con¬ 
tact, m; and, h is the local decrease of the diameter D of a 
particle due to the contact with another particle, m. 

The experimental data in [9, 10] show that the initial 
structure of the powder material corresponds to random 
packing with the parameters Z r = 7.3 and II,. = 0.37 for 
smooth spherical particles (for rough, oblong or irregularly 
shaped particles, especially in highly disperse materials, the 
porosity can reach II = 0.4 or more). We shall assume that 
the unit cells of a material are octahedral (the central particle 
touches four other particles on the top and bottom) with 

II = 0.348 (the average value for two types of octahedral 
packing) but one contact is broken for 70% of the particles. 
Then the coordination number decreases to Z= 7.3, and the 
porosity increases to II = 0.373. The pores form an approxi¬ 
mately cubic structure. 

We shall determine more accurately the structural para¬ 
meters of the initial material, taking account of the nonuni¬ 


formity of the molded material (this is a manifestation of the 
process of clusterization) with a normal distribution of the 
material particles over the coordination number, as follows. 

Let the initial clusters be spherical on average. In accor¬ 
dance with [6] and the experimental data only three-layer 
clusters need be considered. 

In a disperse system, as a result of nonuniform packing 
of the particles, clusters already form as the system is being 
formed (see Fig. 1). During sintering a completely internal 
layer of a cluster consisting of 12 nearest-neighbor particles 
and part of a second, outer layer are formed around particles 
with the maximum coordination number by means of capil¬ 
lary forces. We shall determine the volume fraction of the 
outer layer of the particles in a cluster, which determine the 
cluster and intercluster porosity, as a sum of the correspond¬ 
ing spherical sectors of a three-dimensional system. Let the 
distance / from the center of a cluster to the center of the par¬ 
ticles of the outer layer be the arithmetic mean of the dis¬ 
tance between two rows (layers) in simple cubic packing 2D 
and between two rows in hexagonal [2(2/3)°- 5 ] packing, i.e., 
/= 1.82D. 

Then, from the right triangle OO lt A in Fig. 1 we obtain 
OA = 1.89/7. Next, since the triangles OO n A and OCB are 
similar, we can find the height of the sector (equal to 
D /2 - O n C) of a particle in the outer layer, determining the 
fraction of its volume in the inner part of the cluster. It equals 
approximately 0.43 times the volume of the particle. Using 
the number of neighbors in the first (Z + 2) coordination 
group, we determine the statistical average number N of par¬ 
ticles in the outer layer of the cluster: 

N= 16 + 3Z/4 + Z 2 /8. 

We now use the facts that Z= 12 for the central particle 
of a cluster, Z=(12+Z)/2 for 12 particles in the middle 
layer and the total number of contacts Z £ = 12(12+ Z)/2. 
Then, for the outer layer Z in the cluster part can be taken to 
be the same as in the middle layer Z = (12 + Z)/2. Hence the 
sum of the contacts of the outer layer with the particles of 
the cluster equals Z £ = 0.43 x 0.5(Z + 12)(16 + 3Z/4 + ZV8), 
and the average coordination number Z d for the particles of 
the cluster can be expressed in the following form: 

Z cl = [12 +6(12+ Z) +0.215(12+ Z)(16 + 3Z/4 + 

Z 2 /8)]/[13 + 0.43(16 + 3Z/4 +Z 2 /8)]. (4) 

We shall use the proposed model to analyze the initial 
structure of a green part with point contacts between disperse 
particles with porosity II. 

Using the average porosity and the expressions (1) - (3) 
we shall calculate the average parameters of the structure Z, 
S Q and r. Next, we shall use the normal distribution of the 
particles over the coordination number [6] 

V(Z s = Z) = 48; V{Z s = Z± 1) = 23; V{Z s = Z± 2) = 3 (5) 

for Z < 10. 
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TABLE 1 . Average and Cluster Structure of Unsintered Ceramic Materials 



Average parameters 



Parameters of the cluster structure 


z 

n 

So > 9 m 1 

r, pm 

^rnic Z^mac 

^mic /^mac 

S 0 mic /S 0 mac 

^rnic Z^rnac 

6 

0.476 

0.314 

3.033 

6.280 

0.4598 

0.324 

2.837 

8 

0.373 

0.376 

1.986 

5.597 

8.142 

0.5063 

0.368 

0.296 

0.380 

3.418 

1.937 





7.586 

0.391 

0.365 

2.140 


Given the particle distribution in a cluster, the expres¬ 
sion (5) can be used to calculate the average coordination 
number Z mic of the particles (microstructure) for the inner 
part of a cluster and the average coordination number of the 
particles Z mac (macrostructure) for the intercluster porosity. 

As an example the parameters Z (Z mic and Z mac ), (n mic 
and n mac X S a (S a mic and S 0 mac ) and r (r mic and r mac ) calcu¬ 
lated using the values of n for a system of particles of size 
D = 10 pm (for d = 0) are presented in Table 1. 

It follows from Table 1 that the radii of the capillaries in 
and between clusters differ appreciably. 

During the sintering process the capillary forces in the 
clusters are stronger than in the intercluster space. The 
shrinkage of the cluster formations with the appearance of 
macroparticles in the material is only partial because the 
pores (capillaries) in the intercluster space become larger un¬ 
der the action of tensile stresses. As a result, together with 
shrinkage of the material, sharp increase of the coordination 
number in the clusters (porous particles) Z mic and decrease of 
the microporosity Tl mlc at the start of sintering the average 
size of the intercluster macropores increases. 

We shall now obtain a method for calculating the param¬ 
eters of the porous structure of a material as a function of its 
porosity during sintering, based on the following assump¬ 
tions. 

1. At the same time as the average coordination number 
of the particles in a material increases the sizes of the 
interparticle contacts d/D and, correspondingly, the local de¬ 
formation h/D gradually increase. We assume that h/D in¬ 
creases linearly with decreasing porosity from the sizes of 
the contacts after the formation of a green part h m /D to the 
maximum value h max /D, which for Z = 12 with IT = 0 is ob¬ 
tained using the relation (1) equal to 0.0474. We note that in 
the relation (1) the interparticle contacts are assumed to be 
cylindrical, and as porosity goes to zero the contact surface 
goes to a pentagonal surface. For this reason we shall refine 
Eq. (1) on the basis of the following considerations. 

Two variants are possible as the porosity drops below 

n = 0.26. 

In the first variant a continuous structure comprised of 
elementary dodecahedra, whose unfinished edges form a 
capillary network comprised of connected tetrahedral pores. 
As sintering proceeds the coupled pores transform into iso¬ 
lated pores and can even vanish completely. The result is a 
continuous system of almost regular dodecahedra. The vo¬ 


lume of a spherical particle nD ] / 6 equals the volume of a 
regular dodecahedron with edge a. The volume and surface 
area of the latter are equal to 7.6631a 3 and 20.6457a 2 , re¬ 
spectively. The volume of the dodecahedron can be ex¬ 
pressed as tiZ) 3 /6 = 20.6457a 2 (0.5D - A max )/3. Hence we de¬ 
termine the maximum normal deformation in the contact 
zone as h max /D = 0.0456 at IT = 0. 

The value obtained for h max /D and the value according 
to the expression (1) for Z= 12 and JT = 0 are in the ratio 
0.8718. Then the expression (1) assumes the form 

1 -n = (l -2/Z) 2 /[8( 1.077- l/Z-Z" 116 ) x 

(0.5- 0.8718/i/D) 3 ]. (6) 

In the second variant, as porosity drops below the value 
IT = 0.26 a continuous structure comprised of semiregular 
truncated octahedral is formed. The transformation of the lat¬ 
tice as Z —» 14 into a pore-free structure with n = 0 can be 
represented as compaction of a powder structure owing to 
the interaction of the deformation of the particles with the 
formation of eight primary cylindrical contacts followed by 
six additional secondary contacts. The rectangular through 
channels gradually transform into interconnected tetrahedral 
pores. Subsequent compaction proceeds with the primary 
contacts transforming into hexagonal ones and secondary 
contacts into square ones. The result is a continuous structure 
comprised of particles in the form of semiregular truncated 
octahedra, whose edges form a tetrahedral capillary network 
with Z e = 4. In addition, the capillaries (pores) around square 
contacts are wider than all other capillaries. 

According to the expression (1) the maximum value 
h max /D is 0.04243 for Z= 14 with n = 0. For FI = 0 a dense 
material comprised of semiregular truncated octahedra with 
edge m is formed. They are obtained by cutting off each reg¬ 
ular octahedron with edge a = 3in, whose volume and surface 
area equal 7.6631a 3 and 20.6457a 2 , respectively, six tetrahe¬ 
dral pyramids with edges m. 

The volume of a spherical particle tiD 3 /6 equals the vo¬ 
lume of a semiregular truncated octahedron, which gives 
m = 0.359 D. The heights of the six tetrahedral pyramids cut 
off equal m x 2 0 - 5 and the areas of their bases equal m 2 . 
Hence we obtain that the distance from the center of a poly¬ 
hedron to the hexagonal faces and tetragonal faces equals 
0.4397D and 0.50775 D, respectively. Since the areas of the 
square and hexagonal contacts are in the ratio 1 : 2.5981, we 
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determine the average normal local deformation of the parti¬ 
cles in the contact zone for II = 0 as 

H/D = [8(0.5 - 0.4397) x 2.5981 + 

6(0.5 - 0.50775)]/(8 x 2.5981 + 6) = 0.0451. 

The value of h/D according to the expression (4) for 
Z= 14 and n = 0 equals 0.04245 or 0.9412 times the value 
obtained. Then the refined expression (1) becomes 


4. Let us now examine a unit volume of the material con¬ 
taining a single macroparticle (one cluster). 

The volume and average size of a macroparticle are 
given by 

K mac = ^^ 3 (l + W/6; 

Anac=^(l + ^ ml c)] 1/3 ’ (H) 

where V p mic is the pore volume in a macroparticle, deter¬ 
mined from the expression for the microporosity 


1 -n = (l -2/Z) 2 /[8( 1.077- 1/Z-Z L16 )x 

(0.5 - 0.941A/D) 3 ]. (7) 

Since the correction for Z max =12 (when the particles on 
average are regular dodecahedra) is larger than the correction 
in the expression (7), for low-porosity materials the model 
based on semiregular truncated octahedra is preferable, and 
for this reason we choose it. 

The size of the contacts is set by the porosity of the sys¬ 
tem and the initial value of the average deformation of the 
particles in the contact zone h m /D (or the average initial size 
of the contacts d m /D) after formation: 

h/D = h m /D + 0.0451 - (0.0451 - h m /D) U/U m . (8) 

The average structural parameters Z, S 0 and r are calcu¬ 
lated for the average porosity IT using the expressions (7), 
(8), (2) and (3) for particles of size D. 

2. We assume that before the porosity of ceramic materi¬ 
als starts to decrease during sintering by capillary forces the 
first to form are dense internal regions with spherical (on 
average) particles (microstructure) with average coordination 
number Z mic in the material equal to 

Z mlc = (Z+3)forZ<l\ and Z mic = 14 forZ> 11. (9) 


Hence the micropore volume TI mic in a macroparticle and 
the macroporosity IT mac of the material can be calculated 
from the expressions 

F prmc = ^mic n m,c/( 1 - n m 1 c); n m ac = U ~ V p mic ' ( 12 ) 

5. We assume that h mac /D mac in the macrostructure as¬ 
sumes the maximum value (equal to 0.0451), since the 
macrocontacts are formed by many microparticles in the 
outer layer of a macroparticle. Then the average coordination 
number Z mac can be calculated from the expressions (7). 

Next, the values of S Q and r are determined using the va¬ 
lues of JT mac , Z) mac , h mac /D mac , h\ mac /D mac , d mac /D mac and 
the expressions (2) and (3). 

ALGORITHM FOR CALCULATING 
THE STRUCTURAL PARAMETERS 

1. The average structural parameters are calculated for 
the average porosity II: 

a) the size of the interparticle contacts h/D according to 
the expression (8); 

b) the values of d/D and h 0 /D by iteration according to 
the expressions 


We calculate the microporosity using the expressions (7) 
and (8). The micropore size r mic is determined from the ex¬ 
pression (3) with the values n mic and Z mic for the micro- 
porous space. 

3. To go beyond the cluster sizes given by the expression 
(4) we shall use from it only the initial ratio of the number 
N cl of microparticles forming the microstructure of a cluster 
and the total number of microparticles per unit volume 
containing one average cluster for the initial porosity of the 
material n = 0.373: 

N cl = 13 +0.43(16 + 3Z/4 + Z 2 /8) 
andA j. = 13 + (16 + Z/4 + Z 2 /8), NjN^ = 25.1/43.1 = 0.582. 

The simplest expression for the volume of the material 
corresponding to the internal volume V s mic of a cluster that 
reflects the fractal character of clusterization is 

^mic = ( 1 '-°-7n) 1 ' 8 . (10) 


h = h o -4h;(\.5D-h o )/(3d 2 ) 
and d/D = 2[h 0 /D(l - h 0 /D)]° 5 -, 

c) Z, S 0 and r according to the values of IT, h/D , h 0 /D 
and d/D using the expressions (8), (2) and (3). 

2. The parameters of the microporous structure are calcu¬ 
lated: 

a) the average coordination number Z mic by means of the 
expression (9); 

b) the microporosity Tl mlc by means of the expressions 
(7) and (8); 

c) S Q and r cl by means of the expression (2) and (3) for 
the values of TI mjc , Z mic , h/D , h 0 /D and d/D and the volume 
fraction V s mic of particles in the microporous structure per 
unit volume of the system by means of the expression (10). 

3. The parameters of the macroporous structure are cal¬ 
culated: 

a) the volume V smac and average size D mac of a macro¬ 
particle by means of the expressions (11); 
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TABLE 2. Parameters of the Average, Microporous and Macroporous Structure of Unsintered Ceramic Materials 




Average structure 




Micro- 

and macroporous structures 


n 

z 

h/D 

5 0 ,).im 1 

r. 

According 
to Eq. (6) 

pm 

According 
to Eq. (3) 

^mic /^mac 

^mic Z^mac ^mac > 

So mic /So mac » 

pm - 1 

'ci Aid 

0.37 

7.30 

0 

0.376 

1.981 

1.988 

0.289 

10.3 

36.59 

0.427 

1.354 







0.137 

7.93 


0.141 

1.938 

0.30 

8.43 

0.0088 

0.366 

1.641 

1.760 

0.231 

11.43 

36.23 

0.380 

1.212 







0.104 

8.75 


0.128 

1.619 

0.26 

9.12 

0.0136 

0.352 

1.482 

1.659 

0.197 

12.12 

36.07 

0.347 

1.137 







0.089 

9.19 


0.120 

1.484 

0.20 

10.20 

0.0209 

0.315 

1.272 

1.457 

0.146 

13.2 

35.85 

0.284 

1.030 







0.070 

9.8 


0.104 

1.336 

0.15 

11.19 

0.0268 

0.269 

1.114 

1.261 

0.105 

14 

35.68 

0.220 

0.953 







0.054 

10.36 


0.090 

1.206 

0.10 

12.13 

0.0330 

0.211 

0.949 

1.030 

0.071 

14 

35.44 

0.165 

0.862 







0.033 

11.23 


0.071 

0.919 

0.05 

13.09 

0.0390 

0.138 

0.722 

0.730 

0.036 

14 

35.23 

0.110 

0.663 







0.015 

12.10 


0.050 

0.575 

0 

14.00 

0.0451 

- 

0 

0 

0 

14 

- 

- 

- 







0 

14 





b) the macroporosity Il mac by means of the expression (12); 

c) the coordination number Z mac (by iteration), the spe¬ 
cific surface S a mac and the average pore radius r mac by means 
of the expressions (7), (8), (2) and (3) for the values of n mac , 

^Anac s Anac /^mac ' mac /^inac ttnd t/ mac /D mac . 

The parameters Z (Z mic and Z mac ), n mlc , n mac , D mac , 5 0 
(S 0 mic and S 0 mac ) and r (r mic and r mac ) for a model system of 
particles of size D = 10 pm (d m = 0 for initial porosity 
U r = 0.37) are presented in Table 2. 

It follows from Table 2 that at the start of sintering the 
sizes of the macropores in and between clusters were more 
than 1 /3-rd greater than the pore sizes with uniform 
sintering, after which as the porosity decreases below 15% 
the difference gradually decreases. For porosity below 0.25 
the specific surface in a cluster becomes less that the 
intercluster specific surface. 

CONCLUSIONS 

1. The total and local compaction during sintering of ce¬ 
ramic materials was analyzed and a model for the sintering 
process in such materials and an algorithm for calculating the 
structural parameters for the compaction of these materials in 
the sintering process were developed. 

2. The macropores are more than 1 /4-th larger than the 
pore sizes for uniform sintering. For this reason, the decrease 
of the pore size in a material lags behind the shrinkage of the 
material during heat-treatment, which corresponds to the ex¬ 
perimental data. 
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